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Analytical theory for dark soliton interaction in nonlocal nonlinear materials with 

arbitrary degree of nonlocality 
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We investigate theoretically the interaction of dark solitons in materials with a spatially nonlocal 
nonlinearity. In particular we do this analytically and for arbitrary degree of nonlocality. We 
employ the variational technique to show that nonlocality induces an attractive force in the otherwise 
repulsive soliton interaction. 
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INTRODUCTION 



Spatial optical solitons represent beams, which prop- 
agate in nonlinear media without changing their profile. 
Their existence is a result of an interplay between size- 
determined diffraction and nonlinearity-induced phase 
modulation, which in most cases is produced by the re- 
fractive index modification of the material. Depending on 
the type of nonlinearity, nonlinear media may support ei- 
ther bright or dark solitons [l[ . While bright solitons are 
just finite size beams formed in media with self-focusing 
nonlinearity, dark solitons are more complex objects, as 
they represent an intensity dip in an otherwise constant 
background with nontrivial phase profile I4I . Spatial dark 
solitons have been observed and studied in media with a 
negative or self-defocusing nonlinearity [3|, Ij] . Their tem- 
poral counterparts, which have the form of " dark pulses" , 
i.e., temporal intensity dips on a cw background, can 
exist in optical fibers in the normal dispersion regime 
[5|, M, ■ In recent years the renewed interest in the prop- 
erties of dark solitons stem from experimental advances 
in the physics of matter waves. In particular, the forma- 
tion of dark matter wave solitons have been observed in 
Bose Einstein condensates with repulsive inter-particle 
interaction [7H10|. There has also been report on the 
possibility of dark soliton formation in nonlinear meta- 
materials [10| . Interestingly, temporal dark solitons were 
also shown to be able to induce supercontinuum genera- 
tion in photonic crystal fibers 



II 



The unique property of optical solitons, either bright 
or dark, is their particle- like behavior in interaction [l|. 
However, it is also well known that there is a fundamen- 
tal difference in the interaction of bright and dark soli- 
tons. While bright soltions may attract, repel, or even 
form bound states, depending on their relative phase |12| - 
1J|, dark solitons always repel. This has been confirmed 
in numerous theoretical and experimental works J15l - ll7l | . 
We have shown recently that the nature of dark soli- 
ton interaction can be drastically altered by the spatially 



nonlocal character of nonlinearity 18J, llSf . In nonlocal 
media the nonlinear response of the medium in a par- 
ticular spatial location is determined not only by the 
wave (light) intensity in that position, as in the local me- 
dia, but also by the intensity in a certain neighborhood 
around the point. As a result spatial nonlocality provides 
stabilization of bright solitons |2Cll422l | , and induces their 
attraction, even if they are out-of-phase [231, |2J1 . Nonlo- 
cality has a similar effect on dark solitons. In particular, 
it has been shown both numerically jI8| and experimen- 
tally [19[ that nonlocality induces attraction of otherwise 
repelling dark solitons, leading to the formation of their 
bound states. The physics of soliton attraction in non- 
local nonlinear media can most easily be understood in 
the (linear) regime of strong nonlocality [20| . In the con- 
text of nonlinear optics a strongly nonlocal response of 
the medium leads to the formation a broad (linear) in- 
dex waveguide, which can trap two or more solitons and 
enable the formation of bound states. In the context of 
matter waves such nonlocal (dipolar) interaction leads to 
the formation of a potential well, which again induces 
attraction between solitons. While dark soliton attrac- 
tion has already been observed experimentally [19| the 
theoretical description of this phenomenon in the regime 
of arbitrary degree of nonlocality has been analyzed only 
numerically Il8i or in the special linear regime of strong 
nonlocality j25j . 



In this work we will investigate analytically the interac- 
tion of dark solitons in nonlocal media with an arbitrary 
degree of nonlocality. We will consider a suitable nonlo- 
cal response function and use the variational approach to 
derive analytical formulas for the forces acting between 
two dark solitons. Our results clearly show how nonlo- 
cality induces an attractive force, which depends on the 
degree of nonlocality and counteracts the otherwise in- 
herent repulsive nature of dark soliton interaction. 



THE NONLOCAL MODEL AND THE RESPONSE 
FUNCTION 

In what follows we will be interested in the evolution of 
1+1 dimensional optical beams with a scalar amplitude 
E(x,z) and intensity I{x,z) = \E(x, z)\'^, that depends 
on the transverse x-coordinate and the propagation coor- 
dinate z. Propagation of such beams in materials with a 
nonlocal defocusing nonlinearity can be modeled by the 
following generic nonlocal nonlinear Schrodingcr (NLS) 
equation 



.dE Id^E 



+ 00 



R{x-i)I{^,z)di = Q, (1) 



with the nonlocal response in the form of a convolu- 
tion, where R{x) is the nonlocal response function. In 
what follows we will use the normalization /_ R{x)dx = 
1. Obviously R{x) = 5{x) in a local Kerr medium. 
The actual form of the nonlocal response is determined 
by the details of the physical process responsible for 
the nonlocality. For all diffusion-type nonlinearities 
[26| , orientational-type nonlinearities (like nematic liquid 
crystal) [23|, and for the general quad ratic nonlinearity 
describing parametric interaction |27H30l |. the response 
function is an exponential R{x) = (2cr)~^ exp(— |a;|/cr) 
originating from a Lorentzian in the Fourier domain, with 
(7 > defining the degree of nonlocality. Interestingly, for 
parametric interaction, the response function can also be 
periodic, R{x) oc sin(|a;|/CT) in certain regimes of the pa- 
rameter space. 

To obtain analytically tractable results the strongly 
nonlocal limit of ct — > 00 is often used, in which the 
equation becomes linear 27l|31-33| and the solitons are 
known as accessible solitons [20|. The so-called weakly 
nonlocal limit (cr ^ 1) also presents a simpler model, 
which can be solved exactly for both dark and bright 
solitons [3J|. 

Other types of localized response functions have been 
used to obtain qualitative analytical results that captures 
the physics of the effect of nonlocality, such as a Gaus- 
sian in connection variational calculations 
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The 

generic properties of the different types of response func- 
tions have been studied by Wyller et al. in terms of 
modulational instability and it was shown that in gen- 
eral all types of localized response functions have the 
same generic properties, provided their Fourier transform 
is positive definite j37| . 

Here we combine two approaches. First we use the 
weakly nonlocal model because it allows to study any 
localized response function by a single parameter. This 
allows us to derive the weakly nonlocal form of the inter- 
action potential for any localized response function using 
the variational approach. Then we introduce an arbi- 
trary degree of nonlocality. We do this by assuming a 
box-type localized response function, because this allows 



us to calculate the integrals that appear in the variational 
approach. By comparing the results for arbitrary degree 
of nonlocality and the box-type response to the generic 
results obtained in the weakly nonlocal limit, we prove 
that the results are indeed generic. 



INTERACTION BETWEEN DARK SOLITONS IN 
WEAKLY NONLOCAL MEDIUM 

We begin our analysis by considering first the specific 
weakly nonlocal limit of Eq. ([T|) , in which the width of the 
response function is much smaller than the spatial scale 
of the solitons. Then the intensity of the beam /(^, z) can 
be expanded in a Taylor series with respect to ^ around 
^ = x, and Eq.(l) turns into 



.dE Id^E 



dz 2 dx^ 



EI 



■7 



dx'^ ) 



0, 



(2) 



where 7 = 5 /-_oo R{x)x^dx clearly shows how the re- 
sponse function needs to be localized. It is important 
to note that we have here assumed a symmetric response 
function, which is why it is the second derivative that ap- 
pears as the perturbation term proportional to 7. Asym- 
metric response functions, such as the Raman response 
in optical fibers, could of course easily be used too. How- 
ever, asymmetric response functions do not allow for 
defining a Lagrangian and thus to use the variational ap- 
proach. Thus we consider here only symmetric response 
functions. 

We will investigate the dark solitons using the varia- 
tional (or Lagrangian) approach [38|. It can be shown 
that the Lagrangian density corresponding to Eq.(l2|) is 
of the following form 
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du 
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where we normalized the background intensity of the soli- 
tons to unity and used the following transformation for 
the amplitude of the field E{x,z) = u{x,z)cyi]i(iz). To 
proceed further we must postulate the form of the func- 
tion u{x). It was already shown earlier in studies of local 
dark solitons that the proper ansatz is of the form 



= [B tanh z+ - lA) {B tanh z_ -I- iA) , 



(4) 



where z± = D{x ± xq) and 2xo denotes the separation 
between solitons and A, B satisfy the normalization con- 
dition A^ + B^ = 1. The choice of this particular ansatz 
is dictated by the fact that it represents exact dark soli- 
ton solutions of noninteracting local dark solitons. Sub- 
stituting Eq.(4) into Eq.(3) and considering the case of 



weakly overlapping dark solitons, we obtain the averaged 
Lagrangian L = J_ Cdx in the following form 
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is the Lagran gian for the noninteracting weakly nonlocal 
dark solitons [331 ■ 

From the corresponding Euler-Lagrangian equations 
one finds the following relations for soliton parameters 



^ = 8BDe 
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Assuming well separated and weakly interacting 
(xq-D » 1) almost "black" solitons {A^ k, 0) we can ob- 
tain from Eqs. (7-9) the following equation for the soliton 
coordinate xq 



d'^xo dV{xo) 



(10) 



dz^ dxo 

where we have introduced the "potential" V{xq) as 

F(xo) = yi(a:o) + y2(a:o) (11) 

with 
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FIG. 1: Dark soliton weakly nonlocal interaction potential 
V{xo), given by Eq.(ll). Red dashed line - local case (7 = 0); 
blue solid line - weakly nonlocal regime (7 = 0.05). 
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Therefore, for set values of the soliton parameter B and 
nonlocality 7, the dynamics of soliton interaction is rep- 
resented as a mechanical analogy describing the motion 
of a particle in an external potential. The potential con- 
sists of two contributions. The first one , Vi(a;o), which 
exists even for local nonlinearity, is positive and hence 
is responsible for the naturally occurring dark soliton re- 
pulsion [1^ |40, |4l|. The second contribution, V2(xo), 
provides a nonlocality-mediated attractive force, which 
disappears for 7 = 0. 

The simultaneous presence of competing repulsive and 
nonlocality-induced attractive forces introduces a local 
well in the soliton interaction potential V{xo), as clearly 
demonstrated in Fig. 1 for 7=0.05, which enables the for- 
mation of soliton bound states otherwise not possible in 
the local NLS equation. This result is obtained on the 
specific weakly nonlocal limit, which has the nice advan- 
tage of being generic, in the sense that it is valid for any 
localized and symmetric response function. In the follow- 
ing section we will extend the results to the full regime of 
an arbitrary degree of nonlocality by considering a spe- 
cific response function. However, we connect the general 
results to the generic result of the weakly nonlocal limit 
to demonstrate the generic nature of also the general re- 
sult. 



GENERAL NONLOCAL CASE 

Here we consider the interaction between the dark soh- 
tons in nonlocal media, in which the nonlocal response 
has an arbitrary degree of nonlocality. Then the La- 
grangian density corresponding to Eq.(l) is 



„ i ( ^du du* 
2 \ oz az 



1- 



du 



dx 



1 



+ 00 



- {\u\' -I) j J{x - C) (k(C, z)P - 1) dC- (13) 

In order to make the problem analytically tractable 
we will consider here the particular model of nonlocality 
described by the rectangular nonlocal response function, 



R{x) 



T7- — (T < X < a, 



otherwise. 



(14) 



Physically, this type of nonlocal response means that the 
nonlinear response of the medium in a particular spatial 
location is determined by the equal contributions from 
the light intensity in the neighborhood of this location 
defined by parameter a. This is obviously a simplifica- 
tion, but as we will see later, it leads to a physically 
correct description of the soliton interaction. 

Substituting Eq.(4) into Eq.(12) and integrating over 
transverse coordinate x, we obtain the averaged La- 
grangian in the form 
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Prom the corresponding Euler-Lagrangian equations 
one can derive the evolution equation for the soliton co- 
ordinate, which in the limit of weakly interacting (i.e. 
well separated), almost black solitons (A^ <C 1) takes 
the following form 
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FIG. 2: Dark soliton interaction potential V{xo), given by 
Eq.(16), for a rectangular nonlocal response with an arbitrary 
degree of nonlocality a. Red dashed line - (cr = 1); blue solid 
line - strongly nonlocal regime {a = 8.0). 
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where the effective potential function V{xo) is defined as 
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and parameters B, D and a satisfy the following relation 



1 B^ coth(Dcr) 
3 ~ D2 D^ 



[1 - D^ a^ csch^ (Da)] = 0. (18) 



One can show that in the weakly nonlocal limit, i.e., 
when cr ^ 1 the formula Eq. (16) leads to the potential 
of the form of Eq. (II) with the nonlocality parameter 7 
given by 7 = i J_^ R{x)x'^dx ~ a'^/6. In Fig. 2 we show 
the potential V{xt)) for different values of the nonlocality 
a. It is evident that the generic results of the weakly non- 
local model remain valid also for an arbitrary degree of 
nonlocality, i.e., nonlocality provides an attractive con- 
tribution to the potential, which counteracts the natural 
repulsion of dark solitons thus enabling the formation 



of their bound states. This fact provides evidence that 
our general results for the specific rectangular response 
function are, in fact, generic also for any symmetric and 
localized response function. 

We now confirm our theory by direct numerical simula- 
tions of the nonlocal NLS Eq.(l) with rectangular nonlo- 
cal response function. As initial conditions we used Kerr 
soliton profiles (see Eq.(4)) with A = 0, B = D = 1. 
The representative results are depicted in Fig. 3. These 
contour plots show the dynamics of initially well sepa- 
rated solitons. The separation is chosen in such a way 
that both solitons clearly repel when the nonlinearity is 
local (Fig. 3(a)). It is clear that as the extent of non- 
local response increases both solitons start experiencing 
the attractive force. In fact, in case depicted in Fig. 3(b) 
(cr = 1) the natural repulsion of solitons is almost com- 
pletely compensated for by the nonlocality-mediated at- 
traction leading to the formation of the bound state of 
dark solitons. Interestingly, in this case both solitons arc 
separated by the distance of 2x0=3.8 which corresponds 
to the location of the minimum of the effective potential 
from Fig. 2 for (j=l. For even stronger nonlocality the 
attractive force causes mutual oscillations of solitons tra- 
jectories. The radiation visible in Fig.3(b-d) is a result of 
the fact that the initial wave profiles are not exact dark 
solitons in the nonlocal regime. Hence, the solitons evolve 
and transform as they propagate shading away radiation. 

In Fig. 4 we plot with the solid line the separation 
between solitons corresponding to their bound state as 
calculated from the minimum of the effective potential 
Eq. (|17p . It is evident the separation is nonmonotonic 
function of the degree of nonlocality. This can be ex- 
plained as follows. For small a the nonlocality-mediated 
attractive forces are very weak. Therefore the only way 
to compensate the natural repulsion of the solitons is to 
increase their separation until the latter sufficiently de- 
creases. On the other hand, for large a the nonlocal 
nonlinear potential becomes very broad resulting again 
in an increased separation of solitons. This behavior has 
been confirmed in numerical simulations. To this end, for 
given degree of nonlocality we varied the initial distance 
between the solitons and numerically propagated them 
over distance long enough to establish the formation of 
their bound state. The resulting separation is depicted 
in Fig. 4 by filled squares. Clearly, it follows the trend 
found from variational analysis. On the other hand, the 
numerical data is limited to relatively low degree of non- 
locality because the strong radiation for larger a prevents 
the accurate determination of the bound states. 



CONCLUSION 

We studied analytically the interaction of dark spatial 
solitons in nonlocal medium. We used variational tech- 
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FIG. 3: Numerically simulated interaction of initially parallel 
dark solitons in nonlocal medium with rectangular nonlocal 
response function, for different degree of nonlocality (a) a = 0; 
(b) a — 1, (c) a — 1.5, (d) a = 2.0. Notice the almost par- 
allel propagation of solitons for a = 1 indicating the balance 
between repulsive and attractive forces and, consequently, for- 
mation of the soliton bound state. 



nique to derive the evolution equations for the separation 
between both solitons. We showed that nonlocality pro- 
vides an attractive force between otherwise repelling soli- 
tons. For high degree of nonlocality the attractive force 
is strong enough to enable formation of bound states of 
dark solitons. 
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